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Abstract 


In Fall 2016, ESA presented paper SLS-RFM 16-10 documenting a possible issue with the 
frequency lock-in range specification in Recommendation 2.1.8A of typically 267 to 1067 Hz in 
considerings (b) from considerings (a) for loop bandwidths (2Bz9) in the range of 200 to 800 
Hz with a recommendation of 100 Hz step size for frequency sweeping. The paper calculated 
the lock-in range to be £266 to +1064 rad/s or £42 to +168 Hz. Also, Recommendation 2.1.8B 
has the same issue for considerings (a) and (b), i.e. for 2Bro = 10 Hz, a lock-in range of 13 
Hz was specified and a recommendation of 5 Hz step size for frequency sweeping. ESA also 
provided test results from the Rosetta and Exomars transponders. The results were somewhat 
inconsistent since the tests to verify lock-in and pull-in range did not include acquisition time 
which is vital to the definition of these performance measures. This paper will address these test 
results below. However, we first examine the rationale for Recommendation 2.1.8A/B and its 
consistency with the theory of 2nd order phase lock loop operations. Our approach is to design 
a digital phase locked loop (DPLL) from phase locked loop (PLL) requirements. All analysis 
will be performed with a DPLL. 


1 Introduction 


PLLs originated in 1919 [2]. Since then, the implementation of PLLs into communications hardware 
is commonplace with the DPLL as the preferred platform in modems. Though the origins of 
Recommendations 2.1.8A/B are somewhat obscure, since the input papers are no longer available, 
we assume that the basis for the recommendation originate from classical PLL theory [3]. Therefore, 
we validate these recommendations using a mixture of classical theory and the modern approach 
of DPLL simulations/implementations. We will focus our study on Recommendation 2.1.8A. We 
note that Recommendation 2.1.8B can be extrapolated from these results. 


There are two distinct categories of PLL theory: linear and non-linear analyses. Within each 
category are two more classifications of noise and noiseless analyses. The most well understood 
body of analysis is the linear analysis. However, it can only approximate the acquisition time 


under the noiseless case. On the other hand, noisy non-linear analysis can make good predictions 
for mean-time between cycle slips and the variance of the output jitter, but it cannot predict 
acquisition time. Since Recommendation 2.1.8A is primarily concerned with acquisition, our best 
tool is the noiseless linear theory though the predictions are approximate. However, we note that 
PLLs generally operate with large loop signal-to-noise ratios so effectively they are noiseless. 


This paper has seven sections following the Introduction. Section 2 defines the PLL model and 
how the DPLL was designed. Section 3 discusses the formal linear theory of PLLs and DPLLs. 
Section 4 makes an attempt at non-linear analysis of DPLL. Section 5 defines the PLL performance 
measures from classical linear theory. Section 6 provides the simulation results. Section 7 correlates 
the results with input paper SLS-RFM 16-10. And finally, Section 8 concludes with remarks. 


2 Phase Locked Loop Model 


We first review the design methodology to design a 2nd order DPLL from a 2nd order PLL. 


2.1 Design Methodology 


The approach to DPLL design and simulation is to first start with an analog transfer function /dif- 
ferential equation of the PLL and then transform to discrete time/difference equation of the DPLL. 
This allows the user the ability to use classical control theory parameters, i.e. natural frequency, 
loop bandwidth and damping factor to specify performance characteristics and then simulate/de- 
sign the loop with discrete digital components. Our goal is to design the loop to analyze the lock-in 
range given a particular loop bandwidth based on the natural frequency and damping factor. These 
parameters are used by convention and their definition are found in [1, 3, 2]. The model we will 
use is first defined with s-parameters, i.e. in Laplace transform domain. Figure 1 shows a block 
diagram of this model. The input is a radio frequency sine-wave with frequency of f, and a phase 
offset of 6;. The waveform is passed through a bandpass filter that only delays the signal into the 
detector and does not contribute to the analysis of the loop. The delay effectively adds a value of 7 
to the phase offset. The detector is a multiplier that has a gain defined by Kg. This is followed by 
a control loop filter whose transfer function is F,r(s) which leads to a voltage controlled oscillator 
(VCO) gain followed immediately by the VCO transfer function of ‘. The signal is then fed back 
to the detector. The output of the VCO defines the phase of a cosine function whose frequency is 
the quiescent or expected frequency of the input sine-wave. 


Through trigonometry the multiplier will produce output sine-waves that are a function of the 
phase sum and difference between the input sine-wave and the cosine estimate, i.e. 


sin(27 fot + 6;) x cos(27 fot + 0.) = ; sin(6; — 0.) 4 ; sin(4m fot + 6; + 00) © 5 sin(6, —6,), (1) 


where f, is the frequency of the carrier that is again assumed to be the same as the quiescent 
frequency of the VCO. Note, the frequency of the sum is approximately twice the input frequency 
and will be suppressed since the transfer function of the PLL is a low pass filter. The difference is 
a baseband signal and therefore will determine the operation of the PLL. 


When the input sine wave is locked or acquired, the output of the multiplier is very close to zero 
phase, i.e. 5 sin(6; — 0.) © 5c. This is called the tracking stage and PLL is operating in the linear 
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cos(27 fot + 0) 


Figure 1: Simplified Analog Phase Locked Loop Model 


region and this results in the approximation in (1). In general, during the acquisition mode, the 
PLL is operating in the non-linear region due to the output sine function of the detector. It can 
also be driven into the non-linear region while operating in the linear as a result of noise. For this 
analysis, the linear, noiseless case is only considered. 


2.2 Approach 


The following list of steps define how to design and simulate the loop model using Figure 1: 


1. Define the performance in the analog s-domain (Laplace), i.e. specify the natural frequency 
Wn, and damping factor ¢. Equivalently, the noise loop bandwidth By can be specified since 


2. Set Ka = 1/2 (the scaling value after the multiplication of the input sine function with the 
cosine function from the VCO) and use the active PI loop filter Fy, p(s) = “2+, by defining 


STL 


T2 = 2¢/wy and 7, = 1072. Then Ky = w?27,/Kg. Note steps 1 and 2 are found in [1]. 


3. Use the bilinear transform to convert to discrete-time z-domain, i.e. 


_ 2(z-1) 
°= FGah) a 


where T; is the sampling frequency of the discrete-time PLL!. 


The conversion of F',7(s) to Frr(z) is through the bilinear transform which is now detailed: 


_ Gy + Goz7} 


5 2(z-1) a if = zl 
“Ts (z+1) 


Fr r(z) = Frets) (3) 


where Gy = ifere and Gy = Toor 


Also, the transfer function of the VCO can be defined by combining the VCO gain with the 
integrator and using the bilinear transform: 


‘Note, that because the bilinear transform maps the entire s-plane from [—oo, 0c] on to the unit circle of the 
z-plane [0,27], frequency warping will occur as the transform applies an arctan function to the mapping to avoid 
aliasing. In order to compensate for the distortion, prewarping is needed to ensure that effective equalization of the 
frequency response to largest frequency of interest wo, ie. s = anton [6]. For the active PI DPLL, the 
largest frequencies of interest is at wo = on where 72 = 2¢/wn. Therefore, in our simulations, the largest wo is around 


50 Hz. Then tan (wT; /2) = tan(2%) = 6.25e-4 for T, = 2.5e-5 sec, which is the same as the argument 24 = 6.25e-4. 


Therefore, it is not necessary to use the prewarping transform. 


Fyco(2) = Fvco(s) wo Ie) (4) 


Ts(z+1) 


4. Create difference equations for the loop using (3) and (4). 
5. Simulate the PLL using Matlab. 


3 Formal Analysis Using Linear Theory 


This section starts the analysis using techniques developed in control theory. We begin with the 
classical transfer function analysis. 


3.1 Frequency Transfer Functions 
3.1.1 Continuous Time Transfer Functions 


As a result of restricting the analysis to the linear region, transfer functions can easily be analyzed. 
The open loop transfer function in the s-domain is defined as [3, 1]: 


0,(s) _ KaF ir (s)Ky - Kqwer (stg +1)  w2(s2C/wn +1)  s2Cwy + w2 


= 2 — 2 (5) 


Ho(s) = O.(s) 8 Kqs?7 8 s 


where ©,(s) and @,.(s) are the Laplace transforms of 6, and @ respectively. 
The closed loop transfer function is: 


S82Cwn +w2 


Q.(s)__ Ho(s) eo _ 82un + Wh (6) 


Oi(s) T+ Ho(s) 1 4 nto ~ 5? + 52Cwy + w2 
Ss 


where ©;(s) is the Laplace transform of ;. 


The poles of (6) are at s = —wn(¢+./¢? — 1), found by completing the square of the denominator 
and the zero is at s = — 5. A plot of the poles for various values of ¢ is shown in Figure 2. 
For ¢ < 1, the poles are complex pairs that approach the imaginary axis, i.e. Real{s} = 0. For 
¢ — 0, the response becomes more undamped and oscillatory. These values should be avoided. 
Note that the magnitude of the poles is always w,, for this case, basically saying that the effect on 
the magnitude of the closed loop transfer function (and the error transfer function defined later) is 
stationary at w,. For ¢ = 1, there is a single pole on the real axis at -1. For ¢ > 1, the poles are 


real pairs. For ¢ > 1, the response is purely exponentially damped. 


We can also plot the magnitude frequency response of the closed loop transfer function (6) precisely 
with Matlab for various values of ¢ is shown in Figure 3. 


The phase error transfer function is: 


O-(s) _ 1 _ 1 _ s? 
Oi(s) 1+ Hels) 44 as s* + s2Cwy, + w2’ 
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Figure 3: Closed Loop Transfer Function 


note that since the denominator is the same as in the closed loop transfer function, the pole 
locations for both are identical. We see that there are zeros in numerator existing at location 
s = 0. The response should start at —oco for s = 0 and rise until the locations of the poles at 
Ss = —wy(¢ + \/¢? — 1). In effect, the response should look like a high-pass filter with the knee at 


Wn. A plot of (7) for various ¢ values is shown in Figure 4. 
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Figure 4: Phase Error Transfer Function 


3.1.2. Discrete-Time Transfer Functions 


The denominators of the continuous time transfer functions in (6) and (7), clearly demonstrate the 
second order nature of the system. We can apply the bilinear transform to the loop filter and VCO 


functions in both equations using (2) and get for (6): 


mt Y y26wn + wh 


_ (2-1) AXz—1 AXz—-1 
STD (Fe + fot 2Cwn + wh 


H(z) =]As) 


After some algebra, we get: 


((WnTs)? + 4u%nTsC)z? + 2(WnTs)?z + (wnTs)? — 4unTs¢ 


HAZ) = 


((wnTs)2 + 4wnTs + 4)22 + (2(WnTs)2 — 8)2z + (WnTs)? — 4unTs¢ +4 


For (7), 


2 
(se) 
Ts(z+1) 
H(z) = He (s - — (10) 
eo 2D az—1) \* 4 (2-1) 9 ers 
Ts (z+1) (Sh) T T. (41) Cw PW 
Then, 
Az? —8z4+4 
HAZ) = (nT? 5 5 5 : (11) 
WT)? + 4wyTs¢ + 4) 22 + (2(wn Ts)? — 8)z + (WnTs)? — 4unTs¢ +4 
As in the Laplace domain, H.(z) = 1 — H,(z). 
In Figure 5, we compare the results for H.(f) using (6) by evaluating H,(s) and (8) using 
s=j2rf 
H(z) to show that for the most part, the equations have nearly identical frequency 


z=exp(j27fTs) 

responses. Note that the discrete spectrum will repeat itself every F, = 1/T, Hz due to the 
repetition from the transform z = exp(j27f7T,). An anti-aliasing filter at the digitizer that is not 
shown will suppress the artifacts. However, there is some distortion in the discrete implementation 
before the Nyquist rate, most likely as a result of the frequency warping from the bilinear transform 
[6]. In this case, there is more attenuation of the stop band than predicted by the s-domain. This 
can be seen as beneficial to the overall performance but the transfer characteristics diverge from 
its -20 dB/decade decline. As demonstrated, the s-domain and z-domain responses are on top of 
each other prior to warping as they approach the Nyquist rate. 


We can also plot the phase error transfer functions. In Figure 6, we compare (7) using H-(s) 
s=janf 


and (11) using H.(z) . In this plot, note that the largest suppression occur at frequency 
z=exp(j27fTs) 

below the natural frequency. Above the natural frequency results in phase errors that have little 

or no suppression. Also note that like the closed loop transfer function, there is a zero after the 


Nyquist frequency due to the spectrum repetition. 


3.2 Linear Performance Predictions 


We consider two inputs in the z-domain to (11), i. frequency step of QT: a and ii. frequency 


ramp 01T, 2 ee [2], where Qo is the frequency shift at t = 0 in rad/sec and {; is the constant 


frequency ramp applied at t = 0 in rad/sec? or sec”. 
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Figure 6: Phase Error Transfer Function Magnitude Frequency Response 


3.2.1 Frequency Step Time Response 


For the frequency step: 


6.(z) =0Ts H.(2) 


z 
(2-1)? 
z 4g? — 8244 
(z — 1)? (Wn Ts)? + dun Ts¢ + 4)27 + (2(wnTs)? — 8)z + (Wn Ts)? — 4unTs¢ + 4 
= ADT, z 
~ ((WyTs)? + dary Ts¢ + 4)z? + (2(wy Ts)? — 8)z + (wn Ts)? — 4wnTsC + 4° 


=oTs (12) 


For some convenience, we write the above equation as: 


( eres Zz 
= wn T's +4unTs¢+4 
be(2) a 2 2 (wnTs)?—4 . (@nTs)?—4unTs6+4 3) 

B+ LTT yet donTatt4” + (on Te)2 dan T.CH4 


We can use the z-transform pair: 


245 2d? + b? — 2ab 
az” + bz fe +b ADS scales aa sa) (14) 


22 + 2cz + d? da — 2 


where wg = arccos (—S), ¢ = arctan (<5) and d > 0. For this case, a = 0 then (13) has a 


time response of: 


400T. 
6 (n) = (onTs)?+4u,TitH4 (ae = Aun Ts ne ‘) a cos(wgan ag ?) (15) 
e en eer ae (wpTs)2—4 2 (al)? + 4wy,T.¢ +4 
(wnTs)2+4umTsC+4 ( ta) 
where 
(wnTs)?—4 
ms )2+4wnTsC+4 
Wd = arccos i ~ — — (16) 
, oR 
and ¢ = arctan (—oo) = —7/2. 
Then 
A4Q00Ts 
6(n) = (OnTePanTGTA (‘ent ~ AunTs6 + *) " sinoan) (17) 
e acs (wnTs)2—4 2 (ay Tie)? + 4unTs¢ + 4 
(wnTs)2+4wWnTsC+4 ( aes) 


3.2.1.1 Asymptotic Analysis 


Recall that the final value theorem says that for a discrete function f[k] the limp... f[k] = 
lim,_,1(z — 1)F(z) where F(z) is the z-transform of f[k]. Then (12) becomes, 


zl 
ij ( 1) AN Tz 
= lim (z 
zl ((wnTs)? + 4wnTs¢ + 4)2? + (2(wn Ts)? — 8)z + (wnTs)? — 4unTs¢ + 4 
‘ AQT (27 — z) 
= lim 


21 ((WyTs)? + 4unTs¢ + 4)2? + (2(wn Ts)? — 8)z + (wnTs)? — 4unTs¢ + 4 
—0 
(18) 
Thus the DPLL will have a zero phase error as time approaches infinity for any value of Q9. This 


implies that for a DPLL using an active PI loop filter, the pull-in range, i.e. the frequency range 
for which the DPLL will lock is infinite. 


3.2.2 Frequency Ramp Time Response 


For a frequency ramp, the input is: 0,7) and therefore: 
6.(z) =0,T? = = H(z) 
OT? z(z+1) Az* —8z¢+4 
© 2(z — 1)? ((wn Ts)? + 4unTs¢ + 4)z? + (2(wnTs)? — 8)z + (wnTs)? — 4wnTs¢ + 4 (19) 
0,7? z(z +1) 4 
° 2(z—1) ((wnTs)? + 4wnTs¢ + 4)2? + (2(WnTs)? — 8)2 + (WnTs)? — 4unTs¢ + 4 


20)T? (2? + z) 
(z — 1)(((wnTs)? + 4unTs¢ + 4)2? + (2(WnTs)? — 8)2 + (Wn Ts)? — 4unT5¢ + 4) 


At this point, we can break the equation into partial fractions and find the numerator unknowns 
using residues. Afterwards, we perform an inverse z-transform to each term. Since the development 
is quite tedious, we will not document all steps but to show the end result. 


o% (wnTs)? — 4umTsC +4 1? 
= 2 
6.(n) a | d[n] + un] + A fan + 4a, T C+D) cos(wan + ¢) (20) 
where, 
((nTs)?—4) 
Wq = arccos MGieT a) Wed 4) (21) 
# f (WnTs)2—4wnTsC+4_ |’ 
((wnTs)?+4unTsC+4) 
((wn Ts)? 4am Ts¢—4) (Wn Ts)°—4) (wn Ts)" — dum TsC+4) 
2 w. Ww 
d ee arctan ((wnT's) +4umTs¢+4) ((wn T's) +4wnTs¢+4) (( nT’s) +4 nT's¢+4) (22) 


((wnTs)?—AwnTsC—4)_ | (wnTs)2—4wnTsC-+4 ((wnTs)?—A4) 2] 
((wnTs )? +4 T's¢+4) (wn Ts)? +4onTs ¢+4) ((wnTs )?+4unT's¢+4) 
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and, 


A= 
(Gn Tatum Ta A) |? (wnTs)?—4unTsC+4 | [ en Re yg tern Tacs) |” ((wn Ts)? —4unTs6—4) ((onTs)?—4) (wn Ts)? —4un Ts6+4) 
(en Ts)? +4enTsC+4) | ((wnTs)?+40nTs64+4) " | (wnTs)?+40nTsC+4) (en Ts)? +4enTsC+4) (On Ts)? +40nTsC+4) (Wn Ts)? +40nTs6+4) 
(wn Ts)?—4unTaC+4 (on Ts)?—4) 2 : 
((wn Ts)? +40 T3644) (ois) 
(23) 


Obviously, this is a very unwieldy form of the solution. The insight that can be drawn is that 
the response is in the form of a dampened sinusoid that will rise to final value of oy In general, 
the bilinear transform produces very complex transfer and response functions but it guarantees 
stability in z-domain if there’s stability in the s-domain. We could plot (20) but a simpler approach 
to finding the time response is to use (19) in Matlab and let Matlab perform the inverse. Plots of 
this approach for (20) and (17) respectively, are found in Figures 8 and 7 below. As demonstrated, 
the final values of (20) and (17), i.e. eh and zero respectively, are approached rather quickly, i.e. 


within one lock time for ¢ = 1 and ¢ = 0.707. 


Time response of phase error of 2nd Order PLL with frequency shifts 
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Time response of phase error of 2nd Order PLL with frequency ramps 
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Figure 8: Time Response of 2nd Order DPLL with Frequency Ramp 


3.2.2.1 Asymptotic Analysis 


Using the final value theorem on the frequency ramp response, (19) becomes: 


00 Zz>1 
= lim(z — 1) my | UB (wn Ts)? — dumTs¢ — 4)z + St ((wnTs)? — At Ti +4) 
- Z—1" (WnTe)? + den Tae + Az? + (wns)? — 8)2 + (Wns)? — don Tol +4 
Qy 
we 


(24) 


According to [3], since 0-(k) is actually limited by the sine function, its maximum value is 1 when 
6-(k) = 2/2 = 1.5708 rad. Therefore, 
max Qy = w?. (25) 


Note that (18), (24) and (25) are consistent with the s-domain analysis [3]. 


4 Non-Linear Analysis 


Up to this point, all analysis was performed based on linear methods. Non-linear analysis is 
considerably harder and the results are not insightful. One approach is to use a phase-plane plot 
of the phase error over time. Some work was done on a zero-crossing digital PLL was performed 
in [4]. Though in theory, an argument can be made that there is a direct relationship between 
this DPLL and a zero-crossing DPLL and the results can be interchanged, we will not analyze the 
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zero-crossing DPLL. However, we will introduce the phase-plane technique in this section and make 
use of it in following sections. 


The difference equations that are the result of the DPLL model used can be plotted on a two- 
dimensional Cartesian coordinate system that begins at the initial phase and terminates when the 
DPLL achieves lock. We can either plot the current phase ¢; on the y-axis and the previous phase 
op—1 On the x-axis or plot the difference Ad on the y-axis and the current phase ¢,. Both plots 
contain the same information. Figure 9 shows an example of the technique. Note that all angles are 
normalized to 7 radians. In this case, a frequency step (09/27) of 130 Hz was applied to a DPLL of 
Wry = 120 Hz and ¢ = 0.707. The resulting analysis shows an acquisition time of 0.003575 seconds 
as each plot shows a starting point at around —1.247 radians and an end point at 0 radians. 


4.1 Lock Criteria 


Since a DPLL will vary its estimate when lock is reached, the criteria to determine lock must account 
for this variability. Therefore the criteria that was defined included a measurement accuracy of € 
within zero radians. The most effective lock criteria was the concurrent condition that the current 
and the two previous phase errors all are within ¢ of zero radians. The rationale for this is that 
a second order loop should maintain lock for three consecutive samples. The value for € was 
determined based on the absolute value of 6.(n) (17) and shown in Figure 7 at arrow point based 
on the theoretical lock time defined below (27). The actual € value used deviated slightly within a 
factor of around 2 from (27) since the linear behavior is an approximation. The final values were 
determined by reviewing the actual phase estimation plots, for example reviewing plots such as 
Figure 10 and determining the best values to define lock. Table 1 results are based on € calculated 
values from (27) of -0.048 to -0.063 but the final values used were 0.02 to 0.1. 


The calculated acquisition times were compared to the theoretical lock times and pull-in times. 
The resultant times were very accurate in regards to when the phase error settled to zero. To show 
the efficacy of this approach, Figure 10 shows the phase error analysis used in Figure 9. The input 
phase is shown in black while estimated phase is shown in red. The predicted lock time of 0.0083 
(defined in Section 5.2) is shown in green and the predicted pull-in time (defined in Section 5.4) 
of 0.0014 seconds is shown in blue. The acquisition time of 0.003575 seconds from phase-plane 
analysis is shown in a thicker red line. Note, there is always a phase shift of 7 radians between 
the input phase and the estimated phase due to the delay from the input passband filter. After 
reviewing many simulation runs, the criteria is deemed to be sound and accurate. 


5 Performance Measures Defined 
References [1, 3] provide all of the following performance measures. 


5.1 Lock Range 


The lock range Aw; is “the frequency range within which a PLL locks in one single-beat note 
between reference frequency and output frequency.” 


Aw, % +2¢w,, radians/second. (26) 
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Figure 10: Phase Error Analysis 
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5.2 Lock Time 


The lock time Ty is “the time the PLL needs to get locked when the acquisition process is a (fast) 
lock-in process.” In other words, the lock time is time it takes to lock within the lock range. 


2 
Ty “7 seconds. (27) 


Wn, 


5.3 Pull-in Range 


The pull-in range Awp is “the range which a PLL will always become locked but the process can 
be slow.” For the active PI filter: 


Awp — co radians/second. (28) 


5.4 Pull-in Time 


The pull-in time Tp is “the time the PLL needs to get locked when the acquisition process is a 
(slow) pull-in process.” For the active PI filter: 


120? 


Teo 
P™ 16Cw2 


seconds. (29) 


6 Simulation Results 


6.1 Approach for True Lock Range Determination 


Recommendation 2.1.8A is basically a recommendation of what frequency step to use at the ground 
when sweeping the spacecraft receiver to ensure the steps are within the lock range of the receiver. 
As document SLS-RFM 16-10 points out, the recommendation could be in error in specifying the 
step size in Hz. Our approach is to determine what the true lock range is based on a predicted 
lock time. We will also compare our simulated lock range to the theoretical lock range in equation 
(26). Then we offer conclusions and interpretations of the testing results from document SLS-RFM 
16-10. The simulations are based on a DPLL engine found in [5]. 


To find the true lock range, we find the largest frequency step that will give an average acquisition 
time that is less than or equal to the predicted lock time Tz, (27). We average over eight equally 
spaced incremental initial phase offsets from —7 to 7m. Recall that the acquisition time is found 
using the approach defined in Section 4. 


6.2 Frequency Step Results 


Recommendation 2.1.8A defines the noise loop bandwidth to be in the range of 100 to 400 Hz. For 
a ¢ = 0.707, our goal is to verify that the lock range is defined in (26) or Aw, = £4¢Br/(¢ + i): 
For By = 100 Hz, Aw, = 266.64 radians/second or 42.4370 Hz. The required w, = 2B, /(¢+ i) — 
188.5713 radians/sec or 30.0121 Hz. For By, = 400 Hz, Aw; = 1.0666e+03 radians/second or 


15 


Br 100 Hz 200 Hz 300 Hz 400 Hz 
Natural Frequency 30 Hz 60 Hz 90 Hz 120 Hz 
Theoretical Lock Range +42.4370 Hz | £84.87 Hz | +127.31 Hz | +169.7482 Hz 
Lock Time 0.0333 sec 0.0167 sec | 0.0111 sec 0.0083 sec 
€ 2e-2 4e-2 6e-2 le-1 

Frequency Step 89 Hz 180 Hz 296 Hz 473 Hz 
Pull-in Time 0.0407 sec 0.0208 sec | 0.0167 sec 0.018 sec 
Average Acquisition Times | 0.0333 sec 0.0165 sec | 0.0106 sec 0.0083 sec 

Frequency Step -88.7 Hz -180 Hz -289 Hz -445 Hz 
Pull-in Time 0.0404 sec 0.0208 sec | 0.0159 sec 0.0159 sec 
Average Acquisition Times | 0.0334 sec 0.0166 sec | 0.0105 sec 0.0085 sec 


Table 1: Data Results for ¢ = 0.707 


169.7482 Hz. The required w, = 2Bz/(¢ + iz) = 754.2852 radians/sec or 120.0482 Hz. Lock times 
are Ty, © 2= = 0.0333 seconds for By = 100 Hz and T;, = 0.0083 seconds for B; = 400 Hz. 


Wn 

The objective is use the lock times as a target for the average acquisition times by which we can 
determine the true lock range. Table 1 is a summary of the simulation results. It indicates that the 
theoretical lock range (26) is too narrow. For By; = 100 Hz, a frequency step size of +88 Hz will 
allow the DPLL to lock within the theoretical lock time as opposed to the conservative theoretical 
+42.4370 Hz defined by (26). Also, for By = 400 Hz, a frequency step size of +445 Hz will allow the 
DPLL to lock within the theoretical lock time as opposed to the conservative theoretical +169.7482 
Hz defined by (26). By values for 200 and 300 Hz are provided as well. Each column indicates 
that the lock range prediction (26) is too narrow. The ratios for the By to the true one-sided lock 
range are: 1.1364, 1.1111, 1.0381, and 0.8989 for 100, 200, 300, and 400 respectively. Note that 
the pull-in times are larger than the average acquisition times so their predictions are much less 
useful. 


6.3. Frequency Ramp Results 


For By = 100 Hz and ¢ = 0.707, the result of (25), ie. Q1) = wn = 5.6594e3, does not work. 


Instead a ramp of slightly lower value 0, = pee = 5.3391e3 sec~? will work. A plot of this 
ramp response is shown in Figure 11. Note that the DPLL starts at phase error of 7 radians, i.e. 
the DPLL is in lock since the input passband filter has a delay of a radians. A phase error of 
m +7/2 is reached quickly as predicted by (19) and Figure 8. If we use ¢ = 1, again the value of 


Q) = Oe. = 4.0336e3 sec~? will work (shown in Figure 12) but Q; = wa = 4.07443 sec~? does 
not work. This result also applies to By of 200, 300 and 400 Hz. After, trying various damping 


factors, the conclusion is that (25) is very much an approximation. The DPLL can track a ramp of 
a 


Q) = ears with a <1 for ¢ < 1 anda 2 1 for ¢ > 2. We do not have an explanation for this. We 


2 
speculate that for Q; = 5 to be achieved, the bandwidth of the transfer function has to be large 
enough to accommodate the frequency dynamics. 


Note that it is critical that the DPLL is initially in lock. If the DPLL is not in lock, then our 
results show that the DPLL cannot track even for a = 0.5. 
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Figure 11: Frequency Ramp Response for By = 100 Hz, ¢ = 0.707 and Qy = 5.3391e3 sec~? 
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Figure 12: Frequency Ramp Response for By = 100 Hz, ¢ = 1 and Q; = 4.0336e3 sec~? 


7 Correlating Results with SLS-RFM 16-10 


SLS-RFM 16-10 test results for Section 3.1 are difficult to interpret because there is no acquisition 
time recorded. The definition of lock was a declaration of lock by the receiver. Therefore we cannot 
correlate our results since we don’t know if they achieved a pull-in range or a lock-in range. As 
stated before, the pull-in range for an active-PI loop filter is infinite (28). As a result, we don’t 
understand the test measurements or the results. 


From Table 1, the theoretical lock range for By = 100 and ¢ = 0.707 is +42.4370 Hz and our 
simulation results indicate that true lock range based on achieving the theoretical lock time is +88 
Hz. Given that the SLS-RFM 16-10 tests do not define an acquisition time, no conclusion is drawn 
from comparing results. We find the same conclusion from reviewing section 3.2 of SLS-RFM 
16-10. 


Note, that the sweep rate tests in the annex of SLS-RFM 16-10 to measure w,, will have some 
error given the simulation results in Section 6.3. Also, as noted in the previous section, frequency 
ramping testing must ensure that the PLL start initially in lock or else the PLL cannot maintain 
lock. So the results of calculating w,, from a maximum frequency ramp is problematic as the receiver 
must start in lock. 
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8 Conclusion and Comments 


We are in agreement that Recommendation 2.1.8A is probably in error for considering (b) and 
that recommends (1) should, in principle, be updated. Our simulations indicate that a step size of 
+88 Hz or 176 Hz should be recommended based on the smallest 2B; = 200 Hz. However, since 
the current recommendation of 100 Hz is less than 176 Hz, the recommendation still works but is 
not as efficient. Therefore we do not believe this is not an urgent change. Remember, that our 
simulation results are based on the theoretical lock time of T; = 27/wy. If a different target lock 
time is agreed upon, one that maybe based on a practical time limit, then the results will change. 
We do, however, agree that considering (b) be updated to reflect Hz values immediately. 


As an aside, Gardner [3] page 137, has stated that ¢ < 1 can produce clusters of cycle slips (or 
multiple cycle slips) that can throw the demodulator way out of lock for 2nd order PLL. We have 
seen these cluster slips for Costas loops and DPLLs under simulation and the clusters do reduce in 
likelihood for ¢ > 1. There should be some consideration to recommending ¢ > 1 when missions 
are concerned about cycle slips. If there is interest, we can present our findings in a future input 
paper. 
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